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Abstract 

We present four types of infinitely many exactly solvable Fokker-Planck equations, which are 
related to the newly discovered exceptional orthogonal polynomials. They represent the deformed 
versions of the Rayleigh process and the Jacobi process. 
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I. INTRODUCTION 



The Fokker-Planck equation (FPE) is one of the basic equations widely used for studying 
the effect of fluctuations in macroscopic systems [l]. It has been employed in many areas 
such as astrophysics {2], [J, chemistry {4, 5|, biology ^6 , [7 ] , finance js] and others. Because 



of its broad applicability, it is therefore of great interest to obtain solutions of the FPE for 



various physical situations. Many methods, inc 
ones, have been developed to solve the FPE [l[ 



uding analytical, approximate and numerical 



14j . More recently, the exact solvability of 



the FPE with both time-dependent drift and diffusion coefficients has been studied by the 
similarity method 15]. 

As with any other equation, exactly solvable FPE's are hard to come by. Thus any new 
addition to the stock of exactly solvable FPE's is always welcome. One of the methods of 
solving FPE with time-independent diffusion and drift coefficients is to transform the FPE 
into a time-independent Schrodinger equation, and then solve the eigenvalue problem of the 



latter 



The transformation to the Schrodinger equation of a FPE eliminates the 
first order spatial derivative in the FPE and creates a Hermitian spatial differential oper- 
ator. Owing to this connection, an exactly solvable Schrodinger equation may lead to a 
corresponding exactly solvable FPE. In fact, the well-known solvable FPE describing the 
Ornsterin-Ulenbeck process is related to the Schrodinger equation of the harmonic oscil- 
lator [l|. Another interesting process, called the Rayleigh process, is described by a FPE 



transformable to the Schrodinger equation of the radial oscillator 

Based on this relation, we would like to present infinitely many new exactly solvable 
FPE's which are related to the quantal systems associated with the newly discovered poly- 
nomials, the so-called exceptional orthogonal polynomials ( 19j- 42j). The discoveries of 



these new kinds of polynomials, and the quantal systems related to them, have been among 
the most interesting developments in mathematical physics in recent years. Unlike the clas- 
sical orthogonal polynomials, these new polynomials have the remarkable properties that 
they start with degree t = 1,2..., polynomials instead of a constant, and yet they still form 
complete sets with respect to some positive-definite measure. 

A brief remark on the main developments of these new polynomials is in order. For 



a recent review, see eg. 38|. Two families of such polynomials, namely, the Laguerre- 
and Jacob i-type X\ polynomials, corresponding to I = 1, were first proposed by Gomez- 
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Ullate et al. in Ref. 



19] , within the Sturm-Lioville theory, as solutions of second-order 



eigenvalue equations with rational coefficients. The results in Ref. 



191 ] were reformulated in 



the framework of quantum mechanics and shape-invariant potentials by Quesne et al. 



20 



211 ] . These quantal systems turn out to be rationally extended systems of the traditional ones 



which are related to the classical orthogonal polynomials. The most general Xg exceptiona 
polynomials, valid for all integral i = 1,2,. . ■ were discovered byOdake and Sasaki I 
(the case of £ = 2 was also discussed in Ref. |2l[). Later, in Ref. 23] equivalent but much 



simpler looking forms of the Laguerre- and Jacobi-type Xi polynomials were presented. Such 
forms facilitate an in-depth study of some important properties of the X% polynomials, such 
as the actions of the forward and backward shift operators on the Xi polynomials, Gram- 
Schmidt orthonormalization for the algebraic construction of the Xi polynomials, Rodrigues 
formulas, and the generating functions of these new polynomials. Structure of the zeros of 



24 



40]. Recently, these polynomials were 



21 



28 



29] . An approach which 



the exception polynomials was studied in Ref. 
constructed by means of the Darboux-Crum transformation 
does not rely on Darboux-Crum and shape invariance was reported in 30]. Generalizations 
of these new orthogonal polynomials to multi-indexed cases were discussed in Ref. 34, 



35] . Radial oscillator systems related to the exceptional Laguerre polynomials have been 



considered based on higher order super symmetric quantum mechanics 37 ]. 

So far most of these studies concern mainly with the generation and mathematical struc- 
ture of the new polynomials. Very few physical applications of these polynomials have been 
studied. In this regard, we mention that the role of these new polynomials was considered 



in 25] for the Schrdingier equation with a position-dependent mass, and in Ref. 26 



for 



the Dirac equations coupled minimally and non-minimally with an electromagnetic field. 
Here we show that infinitely many exactly solvable FPE's can be constructed with these 
new polynomials. They represent the generalized, or deformed versions of the Rayleigh and 
Jacobi processes. 

The plan of the paper is as follows. In Sect. II we briefly review the transformation 
between FPE and the Schrodinger equation, and derive the general expression for the prob- 
ability density function of the FPE which is related to the exceptional orthogonal polynomi- 
als. Sect. Ill and IV then discuss the generalized Rayleigh process and the Jacobi process, 
respectively. Sect. V concludes the paper. 
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II. FOKKER-PLANCK AND SCHRODINGER EQUATIONS 



In one dimension, the FPE of the probability density function (PDF) V(x,t) is [H 

V(x,t). (2.1) 



d _.,„. , d 2 



The functions .D^^x) and .D^x) in the FPE are, respectively, the drift and the diffusion 
coefficient (we consider only time-independent case). The drift coefficient represents the 
external force acting on the particle, while the diffusion coefficient accounts for the effect of 
fluctuation. Without loss of generality, in what follows we shall take D^(x) = 1. The drift 
coefficient can be defined by a drift potential U{x) as D^(x) = —U'(x), where the prime 
denotes derivative with respective to x. The current density J(x, t) is given by 

J{x,t) = -^P(x,t)-lp(x,t). (2.2) 



The FPE is closely related to the Schrodinger equation [l|, [16|, |17|]. To show this it is 
convenient to define a function w(x) = —U(x)/2, called the prepotential, so that D^(x) = 
2w'(x). It is the prepotential that occurs naturally in the Schrodinger equation. Substituting 

V(x,t) = e- st e w{x) (j){x) (2.3) 

into the FPE, we find that <fi satisfies the equation Hip = £<p, where 

H = -— + w'{xf + w"{x). 
ox z 

Thus <p satisfies the time-independent Schrodinger equation with Hamiltonian H and eigen- 
value £, and O — exp(w) is the zero mode of H: H<p = 0. 

This implies that a FPE transformable to an exactly solvable Schrodinger equations 
can be exactly solved. One needs only to link the prepotential w(x) in the Schrodinger 
system with the drift coefficient D^(x) = 2w'(x) in the FPE. For example, the shifted 
oscillator potential in quantum mechanics corresponds to the FP equation for the well- 
known Ornstein-Uhlenbeck process 

Specifically, suppose all the eigenfunctions 4> n (n = 0,1,2,.. .) of H with eigenvalues S n 
are solved, then V n (x,t) = e~ £nt (f)o(x)(f) n (x) (n = 0, 1,2, . . .) are solutions of the FPE (12.11) . 
The stationary distribution, corresponding to £q = 0, is Vq(x) = 4>q(x) = exp(2w(x)) (with 
J Vo(x) dx = 1), which is obviously non-negative. Any positive definite initial probability 



density V(x, 0) can be expanded as V(x, 0) = (po(x) ^2 n c n (f> n (x), with constant coefficients 
c n (n = 0, 1, . . .) 

/oo 
n (o;)(0o 1 (^(^O))dar. (2.4) 
-oo 

Note that Co = J™ V(x, 0)dx = 1. At any later time t, the solution of the FPE is 

V{x, t) = (f> (x) ^ c n^n{x) exp(-£ n t). (2.5) 

n 

In this paper we shall be interested in the initial profile given by the 5-function 

V(x,t) = 6(x -x Q ). (2.6) 

This corresponds to the situation where the particle performing the stochastic motion is 
initially located at the point x . From Eqs. (12. 4p and (12.51) we have 

Cn = 0o 1 ( x o)4>n(xo), (2.7) 

and 

V{x,t) = <fi (x)<i) 1 (x )^2<i) n (x )(i) n (x)exp(-£ n t). (2.8) 

n 

We now consider exactly solvable FPEs associated with the newly discovered exceptional 
orthogonal polynomials. Only the singled-indexed exceptional orthogonal polynomials will 
be considered here. There are four basic families of the single-indexed exceptional orthog- 
onal polynomials: two associated with the Laguerre types (termed the LI and L2 types) , 
and the other two with the Jacobi types (termed the Jl and J2 types). Quantal systems 
associated with the Laguerre-type exceptional orthogonal polynomials are related to the 
radial oscillator, while the Jacobi-types are related to the Darboux-Poschl- Teller potential. 



The corresponding FPEs describe the generalized Rayleigh process 18|, and the generalized 
Jacobi process, respectively. 

The eigenfunctions of the corresponding Schrodinger equations are given by 

<J) £ , n (x; A) = N £yn (X)(f) e (x; X)P e , n (v(x); A), (2.9) 

where r](x) is a function of A is cL set of parameters of the systems, A^ >n (A) is a normal- 
ization constant, <pi{x\ A) is the asymptotic factor, and P^ n (a;; A) is the polynomial part of 
the wave function, which is given by the exceptional orthogonal polynomials. When i = 0, 



Po, n {x] X) are simply the ordinary classical orthogonal polynomials of the corresponding 
quantal systems, with P 0i0 (x;X) = constant. When £ > 0, P^ n (x; A) is a polynomial of 
degree £ + n. The PDF (12. 8p for these generalized systems are 

V(x,t) = $(x; X)Pt, {m A)^7oW) 

oo 

x £ N etn (X) 2 P £tn ( V ; X)P iin (rio; X) exp(-£ n (A)t). (2.10) 

n=0 

In the next two sections, we shall consider the Laguerre and Jacobi systems separately. 



III. FPES WITH EXCEPTIONAL X e LAGUERRE POLYNOMIALS 

A. Deformed drift potential and probability density function 

We now consider the exceptional Xg Laguerre polynomials, which appear in the deformed 
radial oscillator potentials. The original radial oscillator potential is generated by the pre- 
potential 

w Q (x;g) = -— + glogx, < x < oo. (3.11) 

The exactly solvable Schrodinger equations associated with the deformed radial oscillator 
potential are defined by the prepotentials (£=1,2,.. .) 

m(x; g) = -^ + (g + £) log* + log U J\ 9+ ^ . (3.12) 
2 6(77; g) 

Here A = g > and rj(x) = x 2 and g) is a deforming function. It turns out there are 



two possible sets of deforming functions £e(r};g), thus giving rise 



to two sets of infinitely 



22 



23l | . These ^ are given 



many exceptional Laguerre polynomials, termed LI and L2 type 
by 

f L^i-n) :L1 
(4 2 \v) :L2 

For both LI and L2 type Laguerre systems, the eigen-energies are £e,n(g) = £n(g+£) = 4n, 
which are independent of g and £. Hence the deformed radial oscillator is iso-spectral to the 
ordinary radial oscillator. The eigenfunctions are given by 

<t>t, n {x',9) = N t>n {g)(j) t ,{x;g)P i n {r];g), <j) £ (x; g) = — — — -, (3.14) 



where the corresponding exceptional Laguerre polynomials Pi >n (r]]g) (£ = 1,2,..., n = 
0, 1,2,.. .) can be expressed as a bilinear form of the classical associated Laguerre polyno- 
mials L"(r]) and the deforming polynomial ^(tj; g), as given in 23j : 



Pi,n{mg) = < 



Uv; g + i)Pn{v; g + £-i)- Uv; g)d v P n (v; g + 1 
(n + g + lr 1 {(g + \Mm 9 + l)P n (w 9 + * + 1) 



1) : LI 



L2. 



(3.15) 



where P n (rj; g) = L n 2 (77). The Xt polynomials P^ n {r\\g) are degree t + n polynomials 
in 77 and start at degree £: Pipiv] g) = ^t{Vi9 + !)■ They are orthogonal with respect to 
certain weight functions, which are deformations of the weight function for the Laguerre 
polynomials (for details, see 23j). From the orthogonality relation of Pi tn (rj;g), we get the 



normalization constants 



23) 



2n\{n+g+l-\) 
(n+g+2l-\)V{n+g+l+\) 

2n\{n+g+\) 
(n + g + t + ±)V(n+ a +l+±) 



: LI 
: L2 



(3.16) 



The drift coefficient and PDF of the FPE generated by the prepotential we(x) in eq.f l3.12p 



are 



x \ = - 2x + 2^ + Ax ( d ^t(vig + 1 ) _ dyUv;g) 



(3.17) 



and 



V(x,t) = (f) 2 £ (x; g)P efi (x 2 ; g)P e ^(x 2 ; g) 

00 

Yl N ln P £,n( x2 > g) P i,n( X l'i 9) exp(-A„t) 



X 



(3.18) 



n=0 



In the limit i — > for g = 1/2, we have &(x; g + 1) — > and g) — > 0, and eqs. (I4.24p 
and (I3.18P become 



-2x + -, 

x 



V(x,t) = 2xe- x2s ^L n (x 2 )L n {xl)e 



Ant 



(3.19) 
(3.20) 



n=0 



where L n (x) are the ordinary Laguerre polynomials. These are just the drift coefficient (with 
the diffusion coefficient being one) and the PDF for the Rayleigh process as given in 18]. 



B. Numerical results 



Below we shall present some numerical results in order to see how the deforming function 
£,i(tj; g) modifies the behavior of the drift- diffusion system. Since the diffusion coefficient 
here is fixed, and the deformation of the system is mainly through the change in the drift 
coefficient, it suffices to see how the drift potential is deformed as I changes. This we plot 
in Fig. 1 for the LI type Laguerre systems. The L2 systems behave similarly, and will not 
be depicted here. Four sets of values of g are shown, namely, g = 0.1,0.5, 1.0 and 5.0. For 
each value of g, we plot Ut(x) with £ = (the undeformed case), 1 and 5. 

Since the drift coefficient is related to the negative slope of Ug(x), one can gain some 
qualitative understanding of how the peak of the probability density moves just from the 
sign of the slope of Ug(x). The reason is as follows. Let us consider the current density 
J(x, t) near the peak, where the slope of the probability density function is very small, 
dP(x,t)/dx ~ 0. From eq. (12. 2p . one has 



As P(x, t) > near the peak, eq. (I3.2ip implies that the peak tends to move to the right (left) 
when it is at a position x such that dUe(x)/dx is negative (positive), until the stationary 
distribution is reached. 

From Fig. 1 one sees that for the values of g considered, the presence of parameter I 
generally makes the slope of Ug(x) more negative near the left wall x = 0. So in the deformed 
systems, the peaks tend to move faster to the right. Also, there can be sign change of the 
slope of the drift potential in certain region near the left wall. In such region, the peak of 
the probability density function will move in different directions for different i. To illustrate 
this, we consider the system with g = 0.5. As mentioned before, the i = case with g = 0.5 
correspond to the well-known Rayleigh process. From Fig. 1 we see that at Xq = 1.2, the 
sign of the slope of U 5 (x) is different from those of U (x) and U\{x). Thus one expects that 
for the initial profile P(x, 0) = 5(x — Xq) with the peak initially located at Xq = 1.2, the peak 
will move to the right for I = 5 system, while for the other two values of £, the peak will 
move to the left. This is indeed the shown in Fig. 2. We note here that the graphs 

shown in Fig. 2 for t = 2.0 are already the stationary distributions of the respective cases. 




peak 



(3.21) 
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IV. FPES WITH EXCEPTIONAL JACOBI POLYNOMIALS 



We now present new exactly solvable FPEs related to the exceptional Jacobi polynomials. 
We shall call these the generalized Jacobi processes. Here we will be brief, as the method of 
construction of these systems is the same as that discussed in the previous section. 

With appropriate redefinition, one can always set the domain of the basic variable to 
be x G [0,7r/2]. The trigonometric Darboux-Poschl- Teller potential has two parameters 
A = (g, h) (g > 0, h > 0) and is defined by the prepotential (£=1,2,.. .) 



7T 



w (x; g, h) = g log sin x + h log cos x, < x < — . 

The corresponding deformation of the system is given by the prepotential 

&{r);g + l,h + l) 



we(x; g, h) — (g + £) log sin x + (h + £) log cos x + log ■ 
where r](x) = cos2x. The deforming function is defined by 



&(V, 9, h) 



(4.22) 



(4.23) 



>(</+*- 



9, h) 



P,'-" ' ' " " ! (//). g>h>0 : Jl 

(-g-t-i, h+t-%). 



P 



'( V ), h>g>0 : J2 



where P^ 1 {x) are the classical Jacobi polynomials. 
The drift coefficients of these systems are 



x) 



2(g + £) cot x-2(h + I) tan x - 2 sin(2x 

The PDF of these systems are again given by eq. ( 12.101) . The basic data needed in the 
formula are 1231 ] : 

£t, n (g, h) = 4n(n + g + h + 2fj, 

<l>t,n{x] 9, h) = (j) e (x; g, h)P t , n (r}] 9, h), M x ', 9, h) 

~2n\ (2n + g + h)Y{n + g + h) 



24) 



(sin x) 9+e (cosx) h+£ 



N n (g,h) 



T{n + g + l)T{n + h + l) 

{n+h+\)(n+g+£-\) 



Nfjg,h)= N 2 n (g + £,h + £) x 



(n+h+l+\)(n+g+2l-\) 

{n+g+\){n+h+i-\) 
{n+g+t+\) (n+h+2£~ |) 



Jl 

J2 



(4.25) 



PeAv,9,h) 



n 



+ h + \y l {{h + g + i,h + i)PX 



(g+e-Lh+e+i) 



(v) 



+ (l+r])Uv,9,h)d v P 



(g+e-Zh+e+i) 



(V)} 



(n + g + ^{(g + \Mm g + hh+ i)pi 9+e+hh+e --*\ v ) 



[ l -v)^(Vi9,h)d v P 



(g+e+ih+i-D 



(V)} 



Jl 



J2 
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As explained in 22|, |23|, the exceptional Jl and J 2 orthogonal polynomials are 'mirror 
images' of each other, reflecting the parity property Pn (—tj) = (—l) n Pn' a \rj) of the 
Jacobi polynomial. 

Below we discuss some numerical results for the J2 systems. In Fig. 3 we plot the drift 
potentials Ue(x) for various sets of parameters. 

Figs. 3(a) and 3(b) show how Ue(x) changes as £ changes at fixed values of g and h. In (a) 
the difference between h and g is small, and the graphs of Ue(x) appear nearly symmetric 
with respect to the mid-point x = 11/ 4, but with the absolute value of the slope becoming 
larger at higher t. So the probability density is pushed faster towards the central region of 
the well at higher £. In (b) h is much larger than g, and one sees that there can be difference 
in the signs of the slopes of the drift potential near the left wall. The effect of the sign 
change is depicted in Fig. 4. 

In Figs. 3(c) and 3(d), we show the change in Ug(x) (at fixed value of £) when either g or 
h is changed while the other is kept constant. It is seen that Ug(x) becomes steeper near the 
left (right) wall as g (h) becomes larger at fixed h (g). This implies the probability density 
function near the left (right) wall is pushed faster towards the right (left) at larger g (h) at 
fixed h (g). 

In Fig. 4 the time evolution of the probably density function is depicted for the three 
cases shown in Fig. 3(b). It has been noted that at some points , the slope of Uf(x) can 
have different signs at different values of £ (at fixed g and h). For instance, at xq = 0.3, 
the slope is positive for I = 0, and are negative for at £ = 10 and 15. If the initial profile 
peaks at xq = 0.3, one expects from the signs of the slopes of Ug(x) that the peak will move 
towards the left wall for the non-deformed system, and Otowards the right for the other two 
deformed cases. Fig. 4 shows that this is indeed the case. 



V. SUMMARY 



We have presented four types of infinitely many exactly solvable FPE's which are the 
generalised, or deformed versions of the Rayleigh process and the Jacobi process. They are 
related to the newly discovered exceptional orthogonal polynomials. These new polynomials 
have the remarkable properties that they start with degree £ = 1,2 ... , polynomials instead 
of a constant, and yet they still form complete sets with respect to some positive-definite 
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measure. In some sense they are deformation of the classical Laguerre and Jacobi poly- 
nomials. The quantal systems that involve these new polynomials are the deformed radial 
oscillator and the trigonometric Darboux-Poschl- Teller potential. Using the transformation 
between FPE and the Schrodinger equation, we obtained the corresponding FPE's, which 
are deformations of the Rayleigh process and the Jacobi process. We have shown numerically 
how the deformation modify the drift coefficient of the FPE, and hence the evolution of the 
probability density function. It gives the possibility to manipulate evolution of stochastic 
processes by suitable choice of a deforming function. 

Acknowledgments 

This work is supported in part by the National Science Council (NSC) of the Republic of 
China under Grants NSC-99-2112-M-032-002-MY3 and NSC-101-2918-I-032-001. The final 
part of the work was done during CLH's visit to the Yukawa Institute for Theoretical Physics 
(YITP), Kyoto University, and he would like to thank the staff and members of YITP for 
the hospitality extended to him, and to R. Sasaki for helpful discussion on exceptional 
orthogonal polynomials. 



11 



[1] H. Risken, The Fokker-Planck Equation (2nd. ed.) (Springer- Verlag, Berlin, 1996). 

[2] S. Nayakshin and F. Melia, Astro-physical Journal Supplement Series 114, 269 (1998). 

[3] G. G. Howes, S. C. Cowley, W. Dorland, G. W. Hammett, E. Quataert and A. A. Schekochihin, 

Astrophysical Journal 651, 590 (2006). 
[4] B. Nowakowski, Phys. Rev. E 53, 2964 (1996). 
[5] Ivan L'Heureux, Phys. Rev. E 51, 2787 (1995). 

[6] M. Schienbein and H. Gruler, Bulletin of Mathematical Biology, Vol. 55, Issue 3, 585 (1993). 

[7] T. D. Frank, P. J. Beek, and R. Friedrich, Phys. Rev. E 68, 021912 (2003). 

[8] R. Friedrich, J. Peinke and Ch. Renner, Phys. Rev. Lett. 84, 5224 (2000). 

[9] R. Blackmorc and B. Shizgal, Phys. Rev. A 31, 1855 (1985). 
[10] A. N. Malakhov and A. L. Pankratov, Physica A 229, 109 (1996). 
[11] T. Blum and A.J. McKane, J. Phys. A 29, 1859 (1996). 
[12] G. Hu, Phys. Rev. A 39, 1286 (1989). 

[13] E. W. Larsen, C. D. Levermore, G. C. Pomraning and J. G. Sanderson, J. of Com. Phy. 61, 
359 (1985). 

[14] H. P. Langtangen, Probabilistic Engineering Mechanics, Vol. 6, 33 (1991). 

[15] W.-T. Lin and C.-L.Ho, Ann. Phys. 327, 386 (2012). 

[16] C.-L. Ho and Y.-M. Dai, Mod. Phys. Lett. B 22, 475 (2008). 

[17] C.-L. Ho and R. Sasaki, Ann. Phys. 323, 883 (2008). 

[18] C.W. Gardiner, Handbook of Stochastic Methods (3rd ed.) (Springer- Verlag, Berlin, 2004). 
[19] D. Gomez-Ullate, N. Kamran and R. Milson, J. Math. Anal. Appl. 359, 352 (2009); 

D. Gomez-Ullate, N. Kamran and R. Milson, J. Approx. Theory 162, 987 (2010). 
[20] C. Quesne, J. Phys. A41, 392001 (2008) ; 

B.Bagchi, C. Quesne and R. Roychoudhury, Pramana J. Phys. 73, 337 (2009). 
[21] C. Quesne, SIGMA 5, 084 (2009). 

[22] S. Odake and R. Sasaki, Phys. Lett. B 679, 414 (2009); 

S. Odake and R. Sasaki, Phys. Lett. B 684, 173 (2009); 

S. Odake and R. Sasaki, J. Math. Phys. 51, 053513 (2010). 
[23] C-L. Ho, S. Odake and R. Sasaki, SIGMA 7, 107 (2011). 



12 



[24] C-L. Ho and R. Sasaki, "Zeros of the exceptional Laguerre and Jacobi polynomials," (to 
appear in ISRN Mathematical Physics, 2012). arXiv: 1102.5669 [math-ph]. 

[25] B. Midya and B. Roy, Phys. Lett. A 373, 4117 (2009). 

[26] C.-L. Ho, Ann. Phys. 326, 797 (2011). 

[27] D. Dutta and P. Roy, J. Math. Phys. 51, 042101 (2010). 

[28] D. Gomez-Ullate, N. Kamran and R. Milson, J. Phys. A 43, 434016 (2010). 

[29] R. Sasaki, S. Tsujimoto and A. Zhedanov, J. Phys. A 43, 315204 (2010) . 

[30] C.-L. Ho, Prog. Theor. Phys. 126, 185 (2011). 

[31] Y. Grandati, Ann. Phys. 326, 2074 (2011). 

[32] S.Odake and R.Sasaki, Phys. Lett. B 682, 130 (2009); 

S.Odake and R.Sasaki, Prog. Theor. Phys. 125, 851 (2011); 
S. Odake and R. Sasaki, J. Phys. A 44, 353001 (2011). 

[33] D. Gomez-Ullate, N. Kamran, and R. Milson, Contemporary Math., 563, 51 (2011). 

[34] D. Gomez-Ullate, N. Kamran and R. Milson, J. Math. Anal. Appl. 387, 410 (2012). 

[35] S. Odake and R. Sasaki, Phys. Lett. B 702, 164 (2011). 

[36] K. Takemura, Heun's equation, generalized hyper geometric function and exceptional Jacobi 
polynomial. arXiv:1106.1543l [math.CA]. 

[37] C. Quesne, Mod. Phys. Lett. A 26, 1843 (2011); 
C. Quesne, Int. J. Mod. Phys. A A26, 5337 (2011). 

[38] C. Quesne, "Exceptional orthogonal polynomials and new exactly solvable potentials in quan- 
tum mechanics", communication at the Symposium Symmetries in Science XV, July 31- 
August 5, 2011, Bregenz, Austria. [arXiv:1111.6467l [math-ph], 

[39] D. Gomez-Ullate, N. Kamran, and R. Milson, A conjecture on Exceptional Orthogonal Poly- 
nomials. arXiv: 1203.6857 [math-ph]. 

[40] D. Gomez-Ullate, F. Marcellan and R. Milson, Asymptotic behaviour of zeros of exceptional 
Jacobi and Laguerre polynomials. larXiv7 l204.2282 [math.CA]. 

[41] B. Midya, Quasi-Hermitian Hamiltonians associated with exceptional orthogonal polynomials. 
larXiv:1205.5860[ math-ph] . 

[42] S. Post, S. Tsujimoto, and L. Vinet, Families of superintegrable Hamiltonians constructed 
from exceptional polynomials. arXiv: 1206.0480 [math-ph]. 



13 





14 



t=0.02 t=0.10 

1.2- 




t=0.5 t=2.0 




1 2 3 4 1 2 3 4 



FIG. 2: Plots of V(x,t) versus x for the LI FPE with different values of t and g = 0.5, xq = 1.2, 
t = (solid curve), 1 (dashed curve) and 5 (dotted curve). Eighty terms in the series in eq. ()3.18|) 
were used. The curves shown in the last sub-figure are the stationary distributions for the respective 
value of I. 
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(c) IrfO, 1=2 (d) g=3, 1=2 




FIG. 3: Plots of Ui{x) versus x for J2 type FPE for : (a) g = 3, h = 4, £ = (solid curve), 1 
(dashed curve) and 15 (dotted curve); (b) g = 1, h = 20, £ = (solid curve), 10 (dashed curve) and 
15 (dotted curve); (c) h = 30, £ = 2, g = 5 (solid curve), 15 (dashed curve) and 25 (dotted curve); 
and (d) g = 3,£ = 2, h = 4 (solid curve), 15 (dashed curve) and 30 (dotted curve) 
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FIG. 4: Plots of V(x, t) versus x for the J2 FPE with different values of t and g = 1, h = 20, xo = 0.3, 
t = (solid curve), 10 (dashed curve) and 15 (dotted curve). Fifty terms in the series in eq. ()2.10|) 
were used. The curves shown in the last sub-figure are the stationary distributions for the respective 
value of £. 
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